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Abstract 

We show that it is possible to describe the effective Pomeron in- 
tercept, determined from the HERA Deep Inelastic Scattering data at 
small values of Bjorken x, using next-to- leading order BFKL evolution 
together with collinear improvements. To obtain a good description 
over the whole range of Q 2 we use a non-Abelian physical renormaliza- 
tion scheme with BLM optimal scale, combined with a parametrization 
of the running coupling in the infrared region. 

1 Introduction & theoretical approach 

The description of Deep Inelastic Scattering (DIS) data for the structure func- 
tion F 2 in different regions of Bjorken x and virtuality of the photon Q 2 is one 
of the classical problems in perturbative QCD. The literature on the subject is 
very large (see, e.g., the reviews in Ref. [I]). In the present letter we are inter- 
ested in regions with low values of x and revisit the theoretical approach to the 
problem using the next-to-leading order (NLO) |2j BFKL [3] equation together 
with collinear improvements. We find that, in order to get a good description 
over the full range of Q 2 , we can use optimal renormalization schemes. In this 
work we focus on indicating which are the most important theoretical aspects 
which drive the bulk of our results. It is possible to introduce subleading 
refinements in our calculation which make our predictions even closer to the 
data and will be presented elsewhere. 
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Let us first review some well-known formulas for DIS in order to set the 
ground for our treatment of the small x resummation. In DIS the cross section 
is written in terms of the structure functions F 2 and F L in the form 

where x and y are the dimensionless Bjorken variables, Q 2 the photon's vir- 
tually and a the electromagnetic constant. More explicitly, for the structure 
functions, in terms of transverse and longitudinal polarizations of the photon, 
we have 

F 2 (x,Q 2 ) = ^-[a T (x,Q 2 ) + a L (x,Q% F L (x,Q 2 ) = ^a L (x,Q 2 ), (2) 



a 



where a^x is the cross-section for the scattering of a transverse (longitudinal) 
polarized virtual photon on the proton. At large center-of-mass energy y/s, 
which corresponds to the small x ~ Q 2 / 's limit, high energy factorization 
makes it possible to write the structure functions Fj, / = 2,L in the form 

where all the integrations take place in the two-dimensional transverse mo- 
menta space with q = The proton ($p) and photon (<3>j) impact factors 
are functions which are dominated by O(Q ) and 0(Q) transverse scales, re- 
spectively. Note that the dependence of $/ on the photon virtuality can be 
calculated in perturbation theory. This is not the case for <3>p whose depen- 
dence on the non-perturbative scale Qo ~ Aqcd can only be modeled. 

If Q 2 was a scale similar to Ql then the gluon Green's function J 7 , which 
corresponds to the solution of the BFKL equation, would be written as 

H ., t , p) = .,' ftf:':K) l i>Y 1 ^ ,4) 



2-rrqp J 2ni J 2-ni \p 2 J \QPj 00 ~ a sXo (7) ' 

with a s = a s N c /ir and Xo(l) = 2^(1) — VKt) — ^(1 — 7) in a leading order 
(LO) approximation, which resums a™ log" s terms to all-orders in the strong 
coupling. ^(7) is the logarithmic derivative of the Euler Gamma function. 
In DIS, however, Q 2 ^> Q\ and this expression should be written in a form 
consistent with the resummation of a s \og{\/ x) contributions: 



2nq 2 J 2m J 2ni \p 2 J \q 2 J u — a s Xo (7 - f) 
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It is well-known that the zeros of the denominator in the integrand generate in 
the limits 7 — > 0, 1 all-orders terms not compatible with DGLAP evolution [1J 
5j. The first of these pieces (0(a 2 )) is removed when the NLO correction to 
the BFKL kernel is taken into account but not the higher order ones, which 
remain and are numerically important. A scheme to eliminate these spurious 
contributions |4), in a nutshell, consists of using a modified BFKL kernel in 
Eq. |4]) where we essentially introduce the change Xo(l) ~~ >" Xo(l + w/2). 

Let us present now in a precise manner our procedure to include the NLO 
corrections and collinear improvements. The NLO expansion of the BFKL 
kernel in terms of poles at 7 = 0, 1 reads 



to 



a s Xo(7 



+ "sXi(7) 



«sXo(7) + a sXi(l) ~ 2 a sXo'(l)xo(l) + 0{a 



a 



s . -2 



7 



+ a£ 



a b 



27= 



+ 



a. 



1 - 

7 27 s 



a: 



2(1 -if 



7 (I-7) 2 2(1 - 7 ) J 



+ 0{a 



s ! ■ 



(6) 



where Xo'(7) = ^'(1 ~~ 7) — ^'{l)- Now, as we have explained before, we resum 
in the Regge region (Q 2 ~ Qq) collinear logarithms by introducing a shift of 
the general form |4l|5] 



CO 



to 



a s (l + Aa s ) 2i[>(l) - V 7 + 77 + 5a s -^1-7+- + ^ 



(7) 



When translated into the DIS limit (Q 2 S> Q 2 ,) this expression is to be replaced 
by 



to = a s (l + Aa s )[2ip(l) -^(7 + Ba s ) + u + Ba s )} 

1 1 2 



a s {l + Aa s )J2 ( 

m=0 ^ 



7 + m + Ba s 1— 7 + m + w-l- Ba s m + 1 



It is possible to get a very good approximation to the solution of this equation 
(certainly within the uncertainty of the resummation scheme) by breaking its 
transcendentality and solving it pole by pole and summing up the different 
solutions. This procedure was proposed in Ref. (5l. In terms of (anti-)collinear 
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poles we obtain 



jr\a s (l + Aa s ) ( 1 



+ Ba s m + 1 

+ - ^7 — 1 - m - Ba s + a/(7 - 1 - m - Ba s ) 2 + 4a s (l + v4a s 
1 1 2 



E 



+ a 



m=0 
2 



a. 



.4 



7 + m 1 — 7 + m m + 1 
+ 
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7 + m 1 — 7 + m (7 + m) 2 (1 — 7 + m) 2 

+ of). 



(1 + m — 7) 3 m + 1 



(9) 



In order to match the NLO poles in Eq. ^ we need to fix A = a and B = —6. 
Keeping the LO and NLO kernels unmodified and introducing only higher 
orders corrections, our collinearly improved BFKL kernel then simply reads 

X(l) = «sXo(7) + «sXi(7) - la 2 s Xo'(l)Xo(l) + Xrg(«s, 7, & )> (10) 



with 



XRG(a s ,7,a ) 6) = ar,(l + aa s ) (^(7) - ip(y - ba s )) 



a 2 



^"(1 — 7) — ba 



^2 



s sin 2 (7r 7 ) ' 2 



2 E I 7 ~ 1 ~ ™ + 

m=0 \ 



2a s (l + aa s 
1 — 7 + m 



+ v% - 1 - m + 6« s ) 2 + 4a s (l + aa s ) . (11) 



For the NLO kernel, 



Xi(7) = 5 Xo(7) 



^-^(7) + 



vl>"(7) + ^ , (l-7)-0(7)-0(l-7) 



7T 2 COS (777) 



4 sin 2 (tt7)(1 -27) 



3 + 



M 2 + 37(1-7) 
+ iV c 3j (3-2 7 )(l + 2 7 ) 



+ ^(3), (12) 



with S = (4 - vr 2 + 5/3o/N c )/12, fa = (^N c - \n f ) and 
0(7) + 0(1 -7) = 



E 

m=0 



1 1 

+ 



7 + m 1 — 7 + m y 



) r ( i+ 2)-* (— j 



(13) 



we obtain for the coefficients 

5 /3o 13 n/ 55 1 (3q rif 11 



(14) 



12 N c 36iV c 3 36' 8iV c 6iV c 3 12' 

As we have already indicated, the non-perturbative proton impact factor has 
to be modeled. We take the functional form 

which introduces three independent free parameters and has a maximum at 
p 2 = 5Qq. Its representation in 7 space reads 

J^ P (p,Ql)(p 2 )-" = 7rCr(5- 7 )(Q2)-7. (16) 

In the present work we choose to keep the treatment of the impact factors 
as simple as possible in order to focus on the behaviour of the gluon Green's 
function. Having this philosophy in mind, we work with the LO photon impact 
factor which can be written in the form (directly inz/ = i(l/2 — 7) space) 

J^i {q, Q 2 ) (^) 7 1 = ««^ 4 E ^ sechM tanh M , (17) 

where Q 2 = (11 + 12z/ 2 )/8 and Q L = v 2 + 1/4. 

So far we have not included in the NLO kernel those terms breaking scale 
invariance, directly linked to the running of the coupling. They appear as a 
differential operator in v space which acts on the impact factors (for a similar 
analysis see Ref. [6]). Exponentiating only the scale invariant LO and NLO 
terms in the kernel, the structure functions can be written as 



Fj(x,Q 2 ) = V / dv x- x \-^ lv ) Cl {v)cp{v) 



1, ■ d , ( ci{v) 



\og{n ) +i— log 



dv \Cp(u 



18) 



where we have gathered different constants in D and \x denotes the renormal 
ization scale at which the QCD coupling is evaluated. Since 

V + V 



Cl (u) = (Q 2 )^ +iu sech(7rz/)tanh(7rz/), (19) 



rrH = i{S- 1 -- i u)(Q 2 r^"\ ,20) 
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we can write our final expression in the form 



Fj(x,Q 2 ) = V I" dvx- x ^ +iy ) Cl {v)c P {v)\ 1 

J — oo I 



log 



Q 2 Q 2 



-i ^7rcoth(7rz/) — 27rtanh (iris) — Mi(y) 



where 



M 2 (u) 



ll + 21z/ 2 + 12z/ 4 
u(l + z/ 2 )(ll + 12z/ 2 ) 



M L (i/) 



" 1 v ( * - - - /?/ 



1 

2 



(21) 



1 - z/ 2 + 4z/ 4 



, (22) 



Although we have included all the ingredients needed to calculate Fl, we leave 
a comparison to experimental data for this observable to future work and focus 
in the following on F 2 . 



2 Running coupling & optimal renormalization 

Although there is some freedom in the treatment of the running of the coupling, 



it is natural to remove the \i dependent logarithm in the second line of Eq. (21 ) 
making the replacement 



«s-«^l°g ^ — > MW, (23) 



, 2 Al og P 

and use this resummed coupling throughout our calculations. We are inter- 
ested in the comparison of our approach with DIS data in the small x region. 
In this letter we focus on the description of the Q 2 dependence of the well- 
known effective intercept X(Q 2 ), which can be obtained from experimental DIS 
data in the region x < 10 -2 through a parametrization of the structure func- 
tion of the form F 2 (x,Q 2 ) = c(Q 2 )x- x(q2 \ The intercept A(Q 2 ) is 0(0.3) at 
large values of Q 2 and 0(0.1) at low values, closer to the confinement region. 
This can be qualitatively interpreted as a smooth transition from hard to soft 
Pomeron exchange. When trying to describe these data with our approach we 
have found that it is rather difficult to get good agreement over the full range 
of 1 GeV 2 < Q 2 < 300 GeV 2 . Somehow it is needed to introduce some new 
ideas related to the infrared region. A recent very interesting possibility is 



6 



that proposed by Kowalski, Lipatov, Ross and Watt in Ref. Alternatively, 
we have found that moving from the MS scheme to renormalization schemes 
inspired by the existence of a possible infrared fixed point significantly helps 
in generating a natural fit for X(Q 2 ), in the sense of having sensible values for 
the two free parameters in our calculation which affect this observable: 5 and 
Qo in the proton impact factor. Here we are guided by having a proton impact 
factor which should be dominated by the infrared region. In the following we 
provide some details of our findings. 

The first evaluation of the BFKL Pomeron intercept in non-Abelian physi- 
cal renormalization schemes using the Brodsky-Lepage-Mackenzie (BLM) opti- 
mal scale setting [8] was performed in Ref. [9] in the context of virtual photon- 
photon scattering. We will use the same procedure in our calculation. The 
pieces of the BFKL kernel at NLO proportional to /3q are isolated and ab- 
sorbed in a new definition of the running coupling in such a way that all 
vacuum polarization effects from the /3q function are resummed, i.e., 



a s (QQo, 7) 



4N r 



A)[log(^) + |xo(7)-| + 2(l + |r)] 



(24) 



where we are using the momentum space (MOM) physical renormalization 
scheme based on a symmetric triple gluon vertex [To] with Y ~ 2.343907 and 
gauge parameter £ = 3 (our results are very weakly dependent on this choice). 
This scheme is more suited to the BFKL context since there are large non- 
Abelian contributions to the kernel. The replacements we need in our kernel 
in order to introduce this new scheme are a s (QQo) —> dt s (QQo) m Eq.(23) 



and Xi(l) ~^ Xi(l) m Eq. (12) together with the corresponding adjustments 
for the coefficients a, b — > a, b which enter Eq. dill) . They read 



Xi{l) = SXoil) 

7i 2 COS (tTj) 

4 sin 2 (tt7)(1 - 27) 
1 

8 



3 *"( 7 ) + *"(!- 7) -0(7) 

2 Uj+ 4 



0(1-7) 



+ 



-CY 
2 y 

13 n/ 

36iV c 3 



6 AT 3 



55 3Y 
"36 H 

11 

12' 



Y 
3" 
3 



16 



1 + 



nf\ 2 + 37(1-7) 
N?) (3 -2 7 )(l + 27) 



3-y 



17Y _ £ 3 
~2 6" 
1 



Xoil), 
17, 



-£ 3 + — Y 
48 s 16 



(25) 
(26) 
(27) 



where S = (4-vr 2 )/12. 
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In order to access regions with Q 2 ~ 1 GeV 2 , we use a simple parametriza- 



tion of the running coupling introduced by Webber in Ref. 11 



a. 



47T 



+ / 



I A 2 



J V A 2 



47T 

Jo 



125 



1 _|_ AIL. 
^ A 2 



1 - VL 
1 A 2 



4+^4 
t -r A 2 



(28) 



At low scales it is consistent with global data of power corrections to pertur- 
bative observables, while for larger values it coincides with the conventional 
perturbative running coupling constant with Landau pole as shown in Fig. [T] 




Figure 1: Left: model for the running coupling with freezing in the infrared (solid 
line) and leading order running coupling with Landau pole (dashed line) for nj = 3 
and A = 0.25 GeV. Right: proton impact factor in momentum space with C = 
l/r(l + 5) and 5, Qq with the values used for the comparison to DIS data. 



The final expression used in the numerical analysis is then given by 



& s (QQo, 7) = &s (QQo, 7) + — / 



(29) 



which replaces Eq. (24) in all expressions. In a future publication we will 



compare the scheme here presented to other physical renormalization schemes. 
For simplicity we have not introduced a complete treatment of quark thresholds 
in the results of this letter, but we have checked that they have a very small 
effect. Let us stress that in our numerical results we do not use any saddle 
point approximation and perform the numerical integrations exactly. 
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3 Comparison to DIS data & scope 



To obtain our theoretical results we have calculated the logarithmic derivative 
Jiogfi/x) usrn § Eq. (21) with the modifications described in Section 2 For the 
comparison with DIS data we chose the values Qo = 0.28 GeV ana 5 = 8.4 
while the dependence on the overall normalization factor C cancels for our 
observable. The QCD running coupling constant is evaluated for rif = 4 and 
A = 0.21 GeV, corresponding to a MS coupling of (M§) = 0.12. The result 
is shown in Fig. [2j The experimental input has been derived from the combined 
analysis performed by HI and ZEUS in Ref. 12 with x < 10 -2 . In the results 
indicated with "Real cuts" we have calculated the effective intercept for F 2 at a 
fixed Q 2 , averaging its values in a sample of x space consistent with the actual 
experimental cuts in x. To generate the continuous line with label "Smooth 
cuts" we have used as boundaries in x space those shown in Fig. |3j which 
correspond to an interpolation of the real experimental boundaries. Note that 
the difference between both approaches is very small. 

We would like to stress the accurate description of the combined HERA 
data in our approach, in particular at very low values of Q 2 . It is noteworthy 
that the values of Qo and S indicate that our proton impact factor (see the 
plot at the right in Fig. [T]) safely lies within the non-perturbative region since 
it has its maximum at ~ 0.81 GeV. In the present letter our intention is to 
emphasize that, in order to reach the low Q 2 region with a collinearly improved 
BFKL equation we needed to call for optimal renormalization and use some 
model with a frozen coupling in the infrared. 

It is possible to improve the quality of our fit by introducing subleading 
contributions such as threshold effects in the running of the coupling, heavy 
quark masses and higher order corrections to the photon impact factor which 



became recently available 13 . We leave these, together with a comparison to 
data not averaged over x, for a more extensive study, which will include an 
investigation of Fl, to be presented elsewhere. 
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Figure 2: Comparison of our prediction with experimental data. 
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Figure 3: Smooth cuts in x used for the effective intercept of F2. 
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